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Abstract 

We analyze long-distance contributions to the K tttt amplitudes relevant 
for the AI = 1/2 selection rule in the framework of the 1/Nc expansion. We 
use a modified prescription for the identification of meson momenta in the chiral 
loop corrections to gain a consistent matching with the short-distance part. Our 
approach involves a separation of non-factorizable and factorizable 1/Nc correc- 
tions. Along these lines we calculate the one-loop contributions from the lowest 
order lagrangian. Our main result is an additional enhancement of the AI = 1/2 
channel amplitude which we find in good agreement with experiment. 



1 Introduction 

Since the first observation of the AI = 1/2 enhancement more than 40 years ago [|I|, 
there have been many attempts to find dynamical mechanisms responsible for it, in 
particular within the standard model. This A J = 1/2 rule was particularly enigmatic 
before the birth of QCD when only the current-current operator Q2 arising from the W 
exchange was considered and, consequently, the ratio ReA{K — > (tttt) j=Q)/ReA{K 
(7r7r)/=2)= ReAo/ ReA2 = R was expected to be about one order of magnitude smaller 
than the experimentally observed value R = 22.2. Now, since the establishment of 
QCD, our understanding of this rule has improved considerably. Within QCD, the 
K TT7T amplitudes are obtained from the effective hamiltonian for AS* = 1 transitions 

Kr = % e ' (1) 

i=l 

involving the Wilson coefficients Ci(/i) which can be calculated for a scale /i ^ 1 GeV 
using perturbative renormalization group techniques, as well as, the local four-quark 

*Prcscntcd at the XXIth School of Theoretical Physics "Recent Progress in Theory and Phe- 
nomenology of Fundamental Interactions" , Ustron, Poland, September 19-24, 1997, published in Acta 
Phys. Pol. B28, 2479 (1997). Two typos in Eq. (17) and Eq. (18) have been corrected. 



operators QiifJ^)- The hadronic matrix elements of these operators are difficult to 
calculate but can be estimated using non-perturbative techniques generally for /i around 
a scale of 1 GeV. For the AI = 1/2 rule only the Zi part of the Wilson coefficient 
Cj is numerically relevant, with Ci(/i) = Zi{fi) — yi{fJ')V^lV-f.^/ {V*gV^^) . The dominant 
operators are 

Ql = 4:SLYdLULlf,UL , Q2 = 4:SL'j''ULUL'Ji,dL , = "8 ^ SlQr QRdL , (2) 

q=u,d,s 

with qji L = |(1 ± 75)g. Major improvements were obtained when it was observed 
that the QCD (and electroweakly) induced effective hamiltonian of Eq.(|l]) can explain 
various large enhancements. These can be of short- distance (SD) nature, like the ffist 
identified octet enhancement in the Q1-Q2 sector dominated by the increase of Z2 
when /i evolves from Mw down to /i ~ 1 GeV. Another important SD enhancement 
was found to arise in the sector of the QCD penguin operators, in particular for zq, 
through the proper inclusion of the treshold effects (and the associated incomplete 
GIM cancellation) 0. Enhancements are also of long-distance (LD) nature like the 
ffist identified LD enhancement of the matrix elements of the QCD penguin operators 
over the matrix elements of Qi and Q2- The latter was ffist conjectured and estimated 
in Ref. using the vacuum insertion method. Due to the non-perturbative character 
of the LD contribution, a large variety of techniques has been proposed to estimate it 
(for some recent publications see Ref. 0). Among the analytical methods, the l/N^ 
approach associated with the chiral effective lagrangian is particularly interesting. 
In this approach, an additional LD enhancement is obtained from inclusion of chiral 
loop effects in the Q1-Q2 sector. The net result of all enhancements mentioned above 
is a value of R in the range of 70-75% of the measured value R = 22.2, suggesting 
that the bulk of the AI =1/2 rule in K ^ tttt decays is now understood.^ One might 
note that the agreement with experiment is not improved by inclusion of the NLO 
renormalization group equations for the Zi^s 0. 

In this proceeding, we reconsider the calculation of ^ tttt amplitudes relevant for 
the AI = 1/2 rule in the 1/Nc approach of Ref. Our main improvement consists 
in the fact that we use a modified matching procedure in order to remove previous 
ambiguities. This will be done treating the factorizable (F) and non-factorizable (NF) 
contributions differently. Only the NF diagrams are matched with the SD contributions 
to cancel the scale dependence of the SD Wilson coefficients (after having, as we will 
see, factorized out the scale dependence of the coefficient 1/m^ in the matrix elements 
of the operator Qe). Factorizable loop contributions which refer uniquely to the strong 
sector of the theory can be calculated in full chiral perturbation theory (xPT), the 
corresponding scale dependence being absorbed in the renormalization of the various 
parameters in the lagrangian. As a result of this procedure, performing a matching 

^In fact the AI = 3/2 channel is usually obtained "sufficiently suppressed" whereas the AI —1/2 
enhancement is obtained only partially. In Ref. |5| for ms(l GeV) — 150 MeV and Aqcd = A'**^ — 300 
MeV the latter channel was reproduced to ~ 70%. 
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with the SD Wilson coefficients, we obtain an additional enhancement of the A/ = 1/2 
channel which we find to be in good agreement with experiment. The A J = 3/2 
channel, however, exhibits a large dependence on the matching scale, resulting from 
the difference of two large numbers, but is found to be equal or smaller than the 
experimental value in such a way that the ratio R = 22.2 is reproduced or even passed 
beyond. 



2 General framework 

To calculate the matrix elements, we start from the chiral effective lagrangian for 
pseudoscalar mesons which involves an expansion in momenta where terms up to 0(p^) 
are included 

/ 



(3) 



with {A) denoting the trace of A and M. = diag(m„, m^, m^). / and r are free param- 
eters related to the pion decay constant (= 93 MeV) and to the quark condensate, 
respectively, with r = — 2(gg)//^. Up to terms of 0{p'^) and to leading order in the 
1/Nc expansion, the lagrangian of Eq.(^ has the most general structure relevant for 
the operators Qi, Q2 and Qq. The degrees of freedom of the complex matrix U are 
identified with the pseudoscalar meson nonet given in a non-linear representation: 



U 



exp -n , 



n = tt'^a. 



{XaXb) = 25, 



n 



where, in terms of the physical states, 

y V2K~ V2K^ 



ab 



V2K+ 



(4) 



(5) 



^a7] J 



with a = cos 9 — v^sin 9 and V2b = sm cos 9. 9 is the rj — rj' mixing angle 

satisfying the relation |^ 



tan26' 



2mlo 



2V2 



m, 



00 



m 



■88 



3a 



1 -1 



ml) 



(6) 



which yields 9 ~ —19°. Note that we treat the singlet as a dynamical degree of freedom. 
Consequently, in Eq.(^) we include the strong anomaly term, with the instanton pa- 
rameter a (~ 0.72 GeV^), which gives a non-vanishing mass of the r/o in the chiral limit 
[rriq = 0) refiecting the explicit breaking of the axial U{1) symmetry. The lagrangian 
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of Eq.(|D is equivalent to the one of Ref. (provided we identify the coefficient 
of Ref. with AL^/ p), except for the fact that we exphcitly include the rjQ. 

A straightforward bosonization yields the chiral representation of the quark currents 
and densities 

f2 

{Jih = q^Ll^qJL = -^^-^{U^^^U),i + ^rL^{^^U^M-M^^^U 

+d^'U^UM^U - U^MU^d^'u) , (7) 

/ ji 

f2 

{Dl\, = {D^j,h = q,Rq,L = -r[^U^ + L,d,U^d>'UU^)^^. (8) 

Using Eqs. (0)-(|l), the operators Qi, Q2 and Qq can be expressed in terms of the meson 
fields. 

The 1/Nc corrections to the matrix elements (Qi)/ are calculated by chiral loop 
diagrams. The loop expansion involves a series in 1//^ ~ which is in direct 

correspondence with the short- distance expansion in terms of as/71 ~ l/A'c. In these 
diagrams we encounter integrals which are regularized by a finite cutoff as it was 
introduced in Ref. p. Due to the truncation to pseudoscalar mesons, the cutoff has to 
be taken at or, preferably, below the (9(lGeV). This restriction is a common feature of 
the phenomenological approaches at hand in which higher resonances are not included. 

To retain the physical amplitudes Aj, which as a matter of principle are scale- 
independent, the long- and the short-distance contributions are evaluated at the cutoff 
scale, i.e., the LD ultraviolet scale is identified with the SD infrared one. Performing 
this identification we must take into account that, within the cutoff regularization, 
there is a dependence on the way we define the momentum variable inside the loop. 

In the standard approach of Refs. |T^, the cutoff is associated to the virtual 

meson, i.e., the integration variable is identified with the meson momentum. Conse- 
quently, as there is no corresponding quantity in the short- distance part, a rigorous 
matching of long- and short-distance contributions is not possible. 

The ambiguity is removed by associating the cutoff to the effective color singlet 



gauge boson as introduced within a study of the K^-Ks mass difference |T2|. This is 
done by identifying in the LD, as well as in the SD part, the momentum variable in 
the loop integration with the momentum flow between the two currents or densities. 
With respect to the standard approach, the momentum of the virtual meson is shifted 
by the external momentum, the former being no longer identical with the integration 
variable, which affects both the ultraviolet, as well as the infrared structure of the 



1/Nc corrections. Note that the matching prescription advocated in Ref. ||T^ was also 



used for current-current operators in the chiral limit in Ref. [13|. The authors showed 



that the coefficient of the quadratic term in the cutoff is increased by a factor 3/2 



relative to the one obtained from the standard matching prescription g 10, 11 1. This 
provides an additional octet enhancement which, however, has to be confirmed by a 
full calculation of the amplitudes relevant for the A/ = 1/2 rule. As we will argue. 
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from this calculation we performed, we confirm this enhancement which is even largely 
increased by the effects beyond the chiral limit. 

Obviously, the modified procedure described above is applicable to the non-factori- 
zable part of the interaction and not to the factorizable ones. The factorizable part, 
however, refers to the strong sector of the theory, and has not to be matched with any 
SD contribution. Consequently, it can be calculated in pure xPT. In this case dimen- 
sional regularization can (and will) be used. Therefore, no momentum prescription 
ambiguities appear. This separation of F and NF contributions was already applied in 
a similar way to investigate the parameter [p^] . 



3 Calculation of the amplitudes and results 

Expanding the lagrangian of Eq.(^, as well as, the currents and densities of Eqs.(|^- 
(H) in terms of pseudoscalars fields, the various F and NF one-loop diagrams to be 
calculated are given in Fig. 1 and Fig. 2 respectively. 






>— ® ® ^ — ^ ® — 
(dq) (qs) ^ (dq) (qs) 

(dq)(qs) (dq)(qs) 



+ (dq)«(qs) 



(dq) (qs) 



Fig. 1. Factorizable diagrams contributing to the matrix elements of the operators Qi- Crossed circles 
represent the currents or densities (with indices (dq) (qs) here specified for Qe)', black circles denote 
strong vertices. The lines represent the pseudoscalar mesons. 



In addition to these diagrams the effect of the wave function renormalization must 
be included, and the values of the various parameters in the lagrangian must be ex- 
pressed in terms of physical quantities. This we did in pure xPT following our general 
procedure for the factorizable contributions. We checked that the scale dependence 
coming from the factorizable diagrams is completely canceled by the scale dependence 
of the various parameters in the tree level expression of the matrix elements. Explicitly, 
we obtain the following expressions: 

^{n^n-mKY = -^{n\'\Q^\Kr = X [l + ^ml) , (9) 

TT 

i{n+n-\Q^\KY = t{n%'>\Q2\KY = , (10) 
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t{n^n-\Q,\KY = ^{n\'\Qe\KY = -A^J-^^)\l , (11) 

with X = V2F^ {m\ — m^) and L5 defined by the relation 

^ = l + ^{ml-ml). (12) 

For the reason of brevity, in Eqs.(|9|)-([TT|) we omit scale-independent (finite) terms re- 
sulting from the one-loop corrections which, nevertheless, are taken into account within 
the numerical analysis. 




Fig. 2. Non-factorizable diagrams with indices here specified for Qq. 

Note that in the definition of L^, in the denominator on the left-hand side of Eq.(|l^), 
we used rather than f. Formally, the difference represents higher order effects. 
Nevertheless, the appearance of / in Eq.(|l2D would induce residual scale dependence 



in Eqs.(|9|)-(|TTD which has no counterpart at the SD level. Therefore, the choice of 
-Ftt is more adequate as it ensures that no scale dependence appears in the complete 
factorizable LD contribution of the matrix elements except the scale dependence of m^. 
To be explicit, in accordance with current conservation, there is no scale dependence 
in Eqs.(^-(|l3); as to the density- density operator of Eq. ([TI|) , we are just left with the 
scale dependence of l/m^. This is to be expected since the evolution of is just inverse 
to the evolution of a single density operator. Consequently, the scale dependence of 
Eq.(|l^) (which is of leading order in A^^^) exactly cancels the corresponding diagonal 
evolution of at SD. This characteristic has already been observed in Ref. ^ . 

The non-factorizable contributions are calculated from the diagrams of Fig. 2, by 
means of introducing a cutoff regulator A Arp. Using the improved matching prescription 
as explained above, we obtain 

i{n'n'\Q,\KY'' = 0, (13) 
z{n^n-\Q,\KY^ = ( " 3 + ^(m^ + 12m^) In A^^) , (14) 
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i{.^.-mKY' = Y^,{l^%F + iml-lml)lnA%,), (16) 
^{.^n'^mKr' = Y6^.rHnA%,, (17) 



Here again we omit finite terms which, however, are taken into account in the numerical 
analysis. The scale-dependent terms in Eqs.(p!3D-(|l8|) have to be matched with the 



Wilson coefficients. To this end, we use the numerical values presented in the leading 



logarithm analyze of Ref. |jT5| 



Note that in the denominators of Eqs.(p!3|)- ([18|) we took the physical value F„ rather 
than /, in the same way as for the factorizable diagrams. Again, the difference repre- 
sents higher order effects. However, the scale dependence of / in Eqs.(|l3D-(|l8]) has no 
counterpart in the SD and will be absorbed at the next order of the chiral expansion. 
Note also that in Ref. only the tree level contribution was taken into account in the 
matrix elements of Qq. The latter are proportional to the L5 coefficient, as the tree level 
contribution of the 0{p^) coming from C^p^) vanishes due to the unitarity of U. From 
the point of the Nc counting it is justified to consider only the tree level contribution for 
Qe, since the Zq coefficient is 0{1/Nc). In Ref. 0, only the one-loop diagrams induced 
by the operators Qi and Q2 were considered. Loops over Qq from CijP'), on the other 
hand, were assumed to be zero as the corresponding tree level contribution from £(p^) 
is zero. However, an explicit calculation of the NF one-loop contributions to the matrix 
elements of coming from C{p^) yields a non-trivial cutoff dependence, as shown in 
Eqs.(p!7|) and (|15), which has to be matched with the SD part. This effect cannot be 
neglected since it corresponds to the leading non-vanishing order of the twofold series 
expansion in l/N^ and p^; i.e., it is of the same order as the tree level term proportional 
to L5 (the former, being of 0{p'^/Nc), is the leading term in the expansion, while 
the latter, being of order 0{p^), is leading in the 1/A^c expansion). The resulting scale 
dependence, however, is only of minor importance within the analysis of the AJ = 1/2 
rule.g 

Our numerical result for the amplitude ReAg is shown in Fig. 3. It shows an addi- 
tional enhancement which renders the amplitude in good agreement with the observed 
value. The new contribution arises from the Qi and Q2 operators. It is due to the 
modified matching prescription in the NF sector (except for approximately one fifth 
of it which is due to the choice of the physical value in Eqs.(|13])-(0) as explained 
above). Our result is remarkably stable with respect to the matching scale. The main 
uncertainty displayed in Fig. 3 originates from the dependence of the Wilson coef- 
ficients on Aqqo (= A*^^)). On the other hand, the isospin 3/2 amplitude shown in 
Fig. 4 is highly unstable. The large uncertainty can be understood from the fact that it 



^ It is actually important within the analysis of the ratio e'/e. Further details on this new calcu- 
lation of the matrix elements of Qq can be found in Ref. |Q6| . 
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is obtained from the difference of two large amplitudes of approximately the same size 
[namely A{K vr+vr") and A{K tt^tt^)]. Consequently, the 3/2 amplitude is not 
well reproduced except that, and this is a crucial point, it comes out to be sufficiently 
suppressed whatever the particular chosen scale is between 500 MeV and 1 GeV. 



MeV 




0J5 - 
0.2 - 

500 550 600 650 700 750 800 850 900 950 1000 

Anf (MeV) 

Fig. 3. Rc^o in units of lO^"* MeV for TOs(1 GeV) = 150 MeV as a function of the matching scale 
AjvF- The experimental value is represented by the dashed line. 

MeV 
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Fig. 4. Rc^2 in units of 10"'^ MeV for ms(l GeV) = 150 MeV. The experimental value is represented 
by the dashed line. 
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In conclusion, it is certainly premature to say that the A/ = 1/2 rule for K — > 
TTTT decays is now understood completely since the l/N^ approach we use is only an 
approximate method. In particular, vector mesons and higher resonances should be 
included in order to take a higher and more secure value for the matching scale. (It 
is probable that the vector mesons play an important role for A2). Nevertheless, we 
believe that the enhancement reported here is a further important indication that the 
1/Nc approach can account for the bulk of the A/ = 1/2 rule for K ^ tctt decays. 

This work has been done in collaboration with G. Kohler, E.A. Paschos and P. 
Soldan. We wish to thank W.A. Bardeen, J. Bijnens, J. P. Fatelo and J.-M. Gerard 
for helpful comments. The author thanks the Deutsche Forschungsgemeinschaft for a 
postdoctoral fellowship in the Graduate Program for Elementary Particle Physics at 
the University of Dortmund. 
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